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THE  USE  OF  THE  DIRECT  MEASUREMENT  OF  SOLAR  RADIATION  FOR  THE 
DETERMINATION  OF  PASQUILL  STABILITY  PARAMETER 


1.  OBJECTIVE 

My  objectives  were  to  compare  two  schemes  for  estimating  the 
amount  of  solar  radiation  received  at  the  earth's  surface  with  solar 
radiation  measurements  and  to  develop  an  algorithm  to  determine  the  Pa s qu i 11 
Stability  Parameter  based  on  direct  measurements  of  solar  radiation. 

2.  BACKGROUND 

Many  organizations,  including  the  U.S.  Army,  have  a  keen  interest 
in  the  rates  of  dispersion  for  vapors,  aerosols,  and  smokes  introduced  to 
the  atmosphere.  These  dispersion  rates  are  needed  to  predict  the  downwind 
concentrations  of  these  substances.  Before  such  estimates  can  be  made, 
however,  atmospheric  turbulence  must  be  understood  and  categorized. 

Atmospheric  turbulence  is  the  key  parameter  in  predicting  concen¬ 
tration  levels.  It  arises  from  two  physical  phenomena,  solar  radiaction 
and  wind.  Solar  radiation  creates  atmospheric  turbulence  because  the 
atmosphere  and  the  earth's  surface  absorb  heat  at  different  rates.  Due  to 
its  physical  properties,  the  earth's  surface  absorbs  heat  much  more  quickly 
than  the  atmosphere.  This  excess  heat  is  transferred  to  the  air  in  contact 
with  the  surface,  which  in  turn  transfers  heat  to  the  layers  of  the  air 
above  it.  This  process  transfers  properties  of  the  lowest  air  layers  to 
those  above  through  mixing,  or  atmospheric  turbulence.  At  night,  the 
surface  cools  faster  and  the  process  reverses  itself. 

The  wind  also  creates  the  mixing  of  air  layers.  As  it  drags 
along  the  (uneven)  surface,  the  wind  creates  small  whirlwinds,  or  eddies, 
that  have  the  same  effect  on  mixing  as  does  solar  radiation.  As  the  wind 
increases,  so  does  the  level  of  mixing.  Therefore,  to  estimate  the  disper¬ 
sion  rates  of  chemicals  released  into  the  atmosphere,  the  combined  effects 
of  solar  radiation  and  wind  on  atmospheric  turbulence  must  be  determined. 

The  Pasquill  Stability  Category  as  revised  by  Turner, 1  is  a  popular 
method  for  characterizing  turbulence  from  the  combined  effects  of  solar 
radiation  and  wind.  An  important  assumption  in  this  procedure  is  that  the 
radiation  reaching  the  earth's  surface  can  be  determined  through  knowledge 
of  solar  altitude,  cloud  cover,  and  cloud  height.  Cloud  height  is  an 
Important  consideration  because  higher  clouds  allow  more  radiation  through 
than  the  lower,  denser  clouds.  Solar  altitude  can  be  determined  precisely 
assuming  a  knowledge  of  latitude,  longitude,  date  and  time.  However,  both 
cloud  cover  and  cloud  height  require  subjective  estimates.  If  direct 
measurement  of  solar  radiation  is  used,  the  need  to  estimate  these  parameters 
would  be  el iminated. 

Routine  meteorological  observations  and  solar  radiation  measurements 
have  been  collected.  The  solar  radiation  measurements  were  compared  with 
subjective  estimates  of  the  same  quantity,  and  a  system  to  determine  the 
Pasquill  Stability  Parameter  based  on  direct  measurement  of  solar  radiation 
is  Included. 
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3. 


SOLAR  RADIATION 


The  energy  in  the  form  of  solar  radiation  received  at  the  top  of  the 
earth's  atmosphere  varies  ver.  little;  the  sun  supplies  a  fairly  constant 
amount.  The  so-called  solar  constant  averages  1340  watts  per  meter  squared, 
or  1.94  langleys  (ly)  per  minute.  The  solar  constant  is  defined  as  the  solar 
radiation  at  normal  incidence  (90°)  outside  the  atmosphere  at  the  mean  solar 
distance  (average  distance  between  the  sun  and  earth).  This  quantity 
fluctuates  as  much  as  ±  1.551,  while  the  varying  distances  between  the  sun 
and  earth  produce  changes  on  the  order  of  ±  3.5%  from  the  mean  value. 2 
However,  solar  radiation  does  not  travel  unimpeded  through  the  atmosphere. 

The  radiant  energy  received  at  the  earth's  surface,  which  determines  atmos¬ 
pheric  turbulence  and  stability,  varies  tremendously  from  place  to  place. 

The  amount  of  radiation  received  at  the  surface  is  dependent  on 
several  factors,  including  latitude,  time  of  day,  cloudiness,  atmospheric 
pollutants,  and  surface  type.  Both  latitude  and  time  of  day  determine  the 
angle  at  which  the  solar  radiation  strikes  the  surface  (radiation  increases 
as  the  angle  Increases). 

Cloudiness  and  atmospheric  pollutants  absorb,  reflect,  and  scatter 
solar  radiation.  Absorption  is  the  process  whereby  solar  radiation  is 
taken  in.  Reflection  involves  the  change  in  direction  of  the  solar  radiation 
after  striking  a  surface,  such  as  a  water  droplet  in  a  cloud.  Scattering 
is  the  mechanism  where  solar  radiation  Is  separated  and  dispersed  in  many 
directions. 

Surface  type  Is  instrumental  In  determining  the  amount  of  radiation 
reflected.  For  Instance,  snow  and  Ice  reflect  much  more  radiation  than  a 
forest-covered  surface.  Table  1  shows  the  approximate  effects  the  atmosphere 
has  on  solar  radiation. 


Table  1.  Solar  Radiation  Lost  to  the  Atmosphere 


Percentage  lost 

Reason 

20 

Absorbed  by  atmosphere  and  clouds 

22 

Reflected  by  clouds 

5 

Scattered  by  atmosphere 

3 

Reflected  by  earth's  surface 

Total  50 

Thus,  only  501  of  the  radiation  reaching  the  top  of  the  atmosphere 
is  received  by  the  earth's  surface.  The  radiation  absorbed  by  a  surface 
unimpeded  by  any  of  these  processes  Is  referred  to  as  direct  radiation. 

The  radiation  reaching  a  surface  as  a  result  of  scattering  or  reflection  from 
another  surface  is  called  diffuse  radiation.  It  Is  interesting  to  note  that 
the  solar  radiation  reaching  the  earth's  surface  can  actually  exceed  the  solar 
constant.  If  the  sun  Is  nearly  directly  overhead  and  the  sky  clear  and 
clean,  very  little  radiation  will  be  absorbed  by  the  atmosphere.  And  if  some 
bright,  fair-weather  cumulus  clouds  are  around,  they  may  add  significant 
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contributions  from  the  radiation  they  reflect.  Thus,  when  diffuse  radiation 
is  considered,  the  earth's  surface  may  actually  receive  more  radiation  than 
the  top  of  the  atmosphere  directly  above. 

4.  PASQUILL-TURNER  STABILITY  CATEGORY 

Tables  2  and  3  are  used  to  determine  the  Pasquill  Stability  Category, 
i.e.,  relative  atmospheric  turbulence.  From  the  existing  meteorological 
conditions,  determine  the  amount  of  the  sky  obscured  by  clouds,  the  approximate 
cloud  height,*  and  the  angle  of  the  sun  above  the  horizon.  Use  Table  2  to 
determine  the  net  radiation  index  (NRI).  Note  that  the  index  ranges  from  0 
to  4  during  the  day  and  -2  to  0  during  the  night.  This  index  reflects 
the  relative  amount  of  solar  radiation  reaching  the  earth's  surface  under 
the  existing  meteorological  conditions.  The  negative  net  radiation  indices 
during  the  night  reflect  the  transfer  of  heat  away  from  the  surface. 

Finally,  the  net  radiation  index  determined  above  is  used  with  the 
wind  speed  to  find  the  Pasquill  Stability  Category,  ranging  from  very 
unstable  (A)  to  very  stable  (F),  in  Table  3.  Tables  2  and  3  summarize  the 
effects  of  solar  radiation  and  wind  on  the  classification  of  atmospheric 
turbulence. 

Some  interesting  aspects  follow  from  a  study  of  Tables  2  and  3. 

For  a  given  NRI,  increasing  the  wind  speed  increases  the  likelihood  that 
the  stability  will  be  D  (neutral).  For  a  given  wind  speed,  increasing  the 
NRI  creates  more  unstable  atmospheric  conditions.  The  greatest  range  of 
stabilities  occurs  at  the  lowest  wind  speeds.  Stabilities  A,  B,  and  C  can 
only  occur  during  the  daytime;  Stabilities  E  and  F  occur  only  at  night;  and 
Stability  D  can  occur  anytime.  In  fact.  Stability  D  is  the  most  common 
atmospheric  condition. 

5.  DATA  COLLECTION 

The  U.S.  Army  Atmospheric  Sciences  Laboratory  (ASL)  has  many  small  , 
specialized  meteorological  teams  located  throughout  the  country,  including 
a  station  near  Phillips  Airfield  at  Aberdeen  Proving  Ground  (APG),  Maryland. 

The  team  routinely  collects  meteorological  data  and  provides  forecasts  for 
various  Army  elements  at  APG. 

Between  the  hours  of  8:00  a.m  and  4:00  p.m,  Monday  through  Friday, 
the  meteorological  team  routinely  collects  such  data  as  temperature,  dew  point, 
pressure,  wind  speed  and  direction,  relative  humidity,  visibility,  precipitation 
cloud  cover,  and  cloud  height.  At  my  request,  the  team  also  measured  solar 
radiation. 


♦Since  cloud  types  are  confined  to  certain  levels  in  the  atmosphere,  cloud 
height  can  be  approximated  reliably  through  a  knowledge  of  cloud  types. 

For  example,  see  Cloud  Charts,  Inc.  Copyright  1958. 


Table  2.  Net  Radiation  Index 


Cloud  cover  (1/8) 

1 _ Day _ , 

Night* 

m-4j 

_ 8 _ 

8 

7-4  3-0 

Cloud  height  (1000  ft) 

i 

>16 1 

16-7 

<7 

<7  >7  | 

Solar  altitude 

■ 

m 

<15° 

fl 

1 

1 

1 

0 

-1  -2 

15-35° 

H 

H 

1 

1 

1 

35-60° 

B9 

Kfl 

2 

1 

2 

>60° 

H 

u 

3 

2 

3 

0 

hour  before  sunset  to  1  hour  after 


sunrise. 


Table  3.  Pasquill  Stability  Category  as  a  Function  of  Net  Radiation 
Index  and  Windspeed 


Wind speed  Net  radiation  index 

m/sec _ 4  3  2  1  0  -1  -2 


<1 

1 

2 

3 

4 

5 

6 

>6 


A  A 
A  B 
A  B 
B  B 
B  C 
C  C 
C  C 
C  D 


B  C 
B  C 
C  D 
C  D 
C  D 
D  D 
D  D 
D  D 


D  F 
D  F 
D  E 
D  E 
D  D 
D  D 
D  D 
D  D 


F 

F 

F 

F 

E 

E 

D 

D 


APG  is  located  in  northeastern  Maryland  on  two  separate  peninsulas 
in  the  northern  section  of  the  Chesepeake  Bay.  The  one-story  meteorological 
building  is  located  approximately  50  meters  from  the  waterfront.  The 
terrain  rises  away  from  the  immediate  area  of  the  building;  i.e. ,  the 
building  and  instruments  nearby  are  located  in  a  localized  low  spot.  In 
addition,  trees  and  other  buildings  almost  completely  surround  the  meteoro¬ 
logical  instrumentation  area.  These  geographical  considerations  no  doubt 
affect  some  of  the  meteorological  data  collected,  particularly  wind  speed 
and  direction. 

Surface  observations  of  the  meteorological  parameters  are  recorded 
between  5  minutes  before  and  the  top  of  each  hour.  All  observations  are 
made  in  accordance  with  the  established  federal  guidelines. ^  The  data  are 
recorded  on  forms  MF1-10A  and  MF1-10B  for  national  quality  control  and 
archiving,  respectively.  For  this  study,  the  meteorological  parameters  of 
most  interest  are  wind  speed,  cloud  cover,  cloud  height,  and  solar  radiation. 
The  solar  elevation  will  also  be  required,  but  these  can  be  determined 
mathematical ly. 

Wind  speed  was  recorded  using  a  Belfort  T420  anemometer  perched 
atop  the  meteorol ogical  building,  approximately  10  meters  above  ground.  Wind 
speeds  are  listed  in  miles  per  hour  and  represent  1-minute  averages.  In 
subsequent  conversations  with  specialists  from  the  meteorological  team,  the 
wind  speed  measurements  are  not  expected  to  be  representati ve  of  the  general 
area.  They  feel  that  the  wind  speed  measurements  recorded  are  lower  than 
should  be  expected  because  the  instrumentation  area  is  lower  than  its 
surrounding  and  the  equipment  was  probably  sheltered  by  the  buildings  and 
trees  and,  in  addition,  the  anemometer  may  not  have  been  properly  calibrated.* 

Cloud  cover  and  cloud  height  are  measured  subjectively.  Cloud 
cover  is  simply  the  amount  of  sky  obscured  by  clouds. 

A  meteorological  specialist  records  cloudiness  through  one  of  four 
descriptors:  clear,  scattered,  broken,  or  overcast.  The  amount  of  cloudiness 
assigned  to  each  category  is  given  in  Table  4. 


Table  4.  Cloudiness  Categories 
Sky  description _ Cloudiness  amount* 


C I  ear 
Scattered 
Broken 
Overcast 


*Scal e  of  0-8 


0 

3 

6 

8 


*Personal  conversation  with  Mr.  Paul  Sisson,  10  August  1984. 


Cloud  height  can  be  determined  quite  accurately  through  a  knowledge  of 
the  cloud  type.  In  general,  cloud  types  are  confined  to  certain  layers  of  the 
atmosphere.  A  cei lometer  is  located  at  Phillips  Airfield,  but  the  instrument 
is  not  often  used.  Also,  on  some  occasions,  balloons  with  known  ascension 
rates  are  released  that  can  determine  cloud  heights.  But,  for  the  most 
part,  cloud  heights  are  estimated  according  to  cloud  types.  Cloud  heights 
are  recorded  in  hundreds  of  feet. 

Solar  radiation  was  measured  employing  a  glass-domed  Eppley  pyra- 
nometer  mounted  1  meter  off  the  ground.  The  device  is  located  approximately 
5  meters  from  the  meteorological  building  and  is  connected  to  it  by  a  ground 
cable.  Solar  radiation  measurements  are  continuously  recorded  inside  the 
building  with  a  pen,  strip  chart,  and  rotating  drum. 

The  pyranometer  measures,  in  langleys  per  minute,  both  direct  and 
diffuse  radiation,  i.e.,  all  the  radiation  being  received  at  that  point. 

The  measurement  of  solar  radiation  was  begun  in  June  1983.  Some 
problems  were  encountered  in  the  initial  weeks  of  data  collection.  The  pen 
used  to  record  the  radiation  on  the  strip  chart  was  often  unreliable, 
including  printing  so  wide  that  the  chart  was  unreadable.  This  was  corrected 
by  the  end  of  June.  Only  data  through  March  1984  judged  to  be  very  good 
quality  was  selected  for  this  study.  The  dates  chosen  are  listed  in  Table  5. 


Table  5.  Dates  Selected  for  Data  Analysis 


Month/ Year 

Dates 

July  1983 

1, 

5,6,7,  8,  18,  19,  20,  21 ,  22 

December  1983 

14, 

15,  16,  19,  20,  21,  22 

February  1984 

14, 

15,  16,  17 

March  1984 

6, 

7,  8,  9 

The  solar  radiation  measurements  received  in  the  form  of  strip 
charts  were  not  directly  usable;  the  data  had  to  be  "reduced"  before  any 
analysis  could  begin.  As  mentioned  previously,  the  strip  chart  provides  a 
continuous  measurement  of  solar  radiation.  Beginning  at  the  top  of  each 
hour,  a  subjective  estimate  of  the  solar  radiation  was  recorded  every  5 
minutes;  i.e.,  an  average  value  for  every  5-minute  increment  was  chosen  to 
represent  the  solar  radiation  strength  over  that  time  interval.  In  addition, 
hourly  averages  were  computed  using  these  5-minute  averages.  For  each 
hourly  average,  the  twelve  5-minute  averages  were  given  equal  weights;  i.e., 
the  12  averages  were  added  and  then  divided  by  12. 

Table  5  represents  a  total  of  25  days  and  214  individual  observations 
of  data.  This  amount  was  judged  large  enough  to  evaluate  the  various  solar 
radiation  estimators  used  in  this  study. 


6. 


SAMPLE  INSOLATION  MEASUREMENTS 


Figure  1  lists  the  solar  radiation  measurements  for  two  dates, 

7  July  and  21  December  1983.  These  dates  were  selected  to  show  the  seasonal 
differences  that  can  be  expected  at  APG.  Table  6  shows  that  7  July  was  a 
mostly  sunny  summer  day,  while  21  December  was  a  mostly  cloudy  winter  day. 


Table  6.  Sky  Conditions  at  Aberdeen  Proving  Ground  on  Select  Dates 

7  July  1983  '  21  December  1983 

Cloud  heights  (ft)/cloud  cover*  Cloud  heights  (ft)/cloud  cover* 


8 

a.m 

Clear 

8  a.m 

12,000  (3) 

25,000  (3) 

9 

a.m 

Clear 

9  a.m 

6,000  (3) 

12,000  (6) 

25,000  (6) 

10 

a.m 

Clear 

10  a.m 

12,000  (3) 

15,000  (6) 

20,000  (6) 

11 

a.m 

Clear 

11  a.m 

8,000  (3) 

12,000  (3) 

15,000  (6) 

12 

a.m 

Clear 

12  a.m 

6,000  (3) 

12,000  (6) 

20,000  (8) 

1 

p.m 

36,000  (3) 

1  p.m 

6,000  (3) 

10,000  (8) 

2 

p.m 

40,000  (3) 

2  p.m 

6,000  (3) 

10,000  (8) 

3 

p.m 

Clear 

3  p.m 

6,000  (3) 

10,000  (8) 

4 

p.m 

Clear 

4  p.m 

6,000  (3) 

8,000  (8) 

*Scale  of  0-8 


The  characteristic  most  evident  in  Figure  1  is  the  difference  in 
total  radiation  received  between  the  2  days.  Some  of  this  difference  is 
explained  by  the  cloudiness  on  21  December,  but  most  may  be  explained  by 
two  factors.  APG  receives  radiant  energy  over  a  longer  time  period  (approxi¬ 
mately  5  hours)  in  July  than  it  does  in  December  and,  secondly,  APG  receives 
more  direct  radiation  in  summer  than  it  does  in  winter.  In  this  case,  a 
surface  at  APG  on  21  December  1983  would  have  received  only  13*  of  the 
total  radiation  it  would  have  received  on  7  July  1983. 

The  peak  value  of  1.39  langleys  per  minute  on  7  July  1983  occurred 
at  1:15  p.m,  while  the  corresponding  value  for  21  December  1983  of  0.48 
langley  per  minute  occurred  at  11:30  a.m. 

7.  DATA  SUMMARY  -  NET  RADIATION 

Solar  radiation  measurements  have  been  compared  to  various  subjecti 
estimates  of  that  quantity.  For  purposes  of  this  study  the  radiation 
measurements  were  reduced  in  three  forms:  an  "instantaneous"  measurement 
at  the  top  of  each  hour,  and  the  average  hourly  values  previous  to  and 
following  this  time.  The  instantaneous  measurement  actually  represents  a 
subjective  estimate  of  the  average  5-minute  solar  radiation. 

Table  7  lists  the  geometric  means  for  these  measurements  of  solar 
radiation  for  each  net  radiation  category.  The  total  number  of  occurrences 
do  not  total  214  since  some  occurrences  of  NRI  equaling  -1  and  -2  existed. 
Recall  that  radiation  data  was  collected  near  sunrise  in  winter  months. 
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This  time  period  is  defined  as  ’'night”  by  the  Pasquill-Turner  (P-T)  System, 
and  the  stability  category  was  often  E  or  F.  These  10  data  points  have 
been  eliminated.  Only  Pasquill  Stability  Categories  A  through  D  were 
analyzed. 


Table  7. 

Mean  Values* 
0  through  4 

of  Solar 

Radiation  for  Net 

ladiation  Categori 

es 

8 

£ 

Number  of 

Hourly 

Number  of 

Hourly 

& 

NRI 

occurrences 

average 

Instantaneous 

occurrences 

delay 

| 

0 

45 

.16 

.15 

38 

.18 

% 

1 

43 

.30 

.33 

39 

.33 

;S; 

2 

47 

.49 

.52 

45 

.52 

.v 

3 

42 

.83 

.85 

33 

.90 

4 

27 

1.08 

1.08 

27 

1.14 

fl 

204 

182 

o 
,  >  1 

♦Langleys  per  minute 

Table  7  shows  that  the  P-T  net  radiation  index  (NRI)  is  intrinsically 
correct.  As  the  NRI  increases,  so  do  the  solar  radiation  measurements. 

Note  that  all  three  measurements  of  solar  radiation  are  comparable  with  one 
another. 

The  geometric  means  given  in  Table  7  list  the  approximate  value  for 
each  NRI  category,  but  it  presents  nothing  in  relation  to  the  distribution 
of  the  solar  radiation  measurements  about  the  means.  Tables  8  and  9  list 
the  joint  occurrence  of  the  NRI  category  and  solar  radiation. 

Table  8.  Joint  Occurence  of  NRI  and  Instantaneous  Radiation  (204  Observations) 

Instantaneous  insolation* 


NRI 

<.2 

.20-. 39 

.40-. 69 

.70-1.09 

>  1.1 

Total 

0 

16.2% 

4.4 

0 

1.5 

0 

22.1 

1 

7.8 

9.3 

2.5 

1.0 

0.5 

21.1 

2 

1.9 

5.9 

11.3 

3.4 

0.5 

23.0 

3 

0 

2.5 

2.9 

10.8 

4.4 

20.6 

4 

0 

0.5 

1.5 

2.9 

8.3 

13.2 

25.9 

22.6 

18.2 

19.6 

13.7 

1U0.0 

Table  9.  Joint  Occurence  of  NR1  and  Hourly  Delay  Radiation  (182  Observations) 


m 

.20-. 49 

>  1.2 

Total 

0 

14.8* 

4.9 

0.6 

0.6 

0 

20.9 

1 

7.7 

9.9 

3.3 

0.6 

0 

21.5 

2 

2.7 

7.2 

13.2 

1.6 

0 

24.7 

3 

0 

1.1 

4.9 

10.4 

1.6 

18.0 

4 

0 

0.6 

1.1 

5.5 

7.7 

14.9 

757? 

23.7 

23.1 

18.7 

9.3 

100.0 

*Langleys  per  minute 


Table  8  is  a  comparison  of  the  instantaneous  value  of  the  solar 
radiation  at  the  same  moment  the  meteorological  observations  were  being  made 
to  determine  NRI.  Table  9  shows  how  well  the  NRI  value  predicts  the  average 
solar  radiation  over  the  next  hour.  However,  the  results  show  that  the 
Pasquill-Turner  (P-T)  system  did  not  fare  well  in  this  comparison.  For 
instance,  the  sum  of  the  observations  along  the  diagonal  for  Tables  8  and  9 
total  55.9  and  56. OX,  or  barely  more  than  half  the  time.  If  there  were  a 
perfect  relationship  between  these  two  systems  of  measurement,  all  the 
observations  would  be  found  along  the  diagonal. 

8.  DATA  COMPARISON  -  INSOLATION  PARAMETER  0 

An  interesting  approach  to  estimating  solar  strength  that  relies 
somewhat  on  subjective  measurements,  but  not  as  much  as  the  P-T  system,  is 
described  below.  Eauation  1  describes  the  flux  density  of  solar  energy  on  a 
horizontal  surface. 4 

Solar  radiation  strength  =  K  (1  -  AN)  sin  a  (1) 


where , 

K  =  proportionality  factor  dependent  on  the  solar  constant  and  atmospheric 
transmission  (the  amount  of  radiation  reaching  the  earth's  surface 
allowing  for  losses  due  to  atmospheric  scattering  and  absorption) 

A  =  average  cloud  albedo  (the  percentage  of  incident  radiation  reflected) 

N  =  percentage  of  sky  covered  by  clouds 

a  =  solar  elevation  angle  (above  the  horizon) 

Equation  1  applies  to  an  area  sufficiently  large  that  the  consequences  of  the 
irregular  distribution  of  clouds  are  minimized.  If  K  were  assumed  to  be 
constant  and  an  appropriate  value  for  A  chosen,  such  as  0.5,  an  insolation 
parameter  t  can  be  defined: 

t  =  (1  -  0.5  N)  sin  a  (2) 

where , 


*  ranges  between  0  and  1. 


The  solar  elevation  at  any  point  on  the  earth's  surface  can  be 
estimated  very  accurately  using  equation  3. 

sin  a  =  sin  e  sin  6  +  cos  6  cos  6  cos  h  (3) 

where, 

e  =  latitude  at  the  point  of  observation 

6  =  declination  of  the  sun  (angular  distance  north  or  south  of  the  celestial 
equator) 

h  =  hour  angle  of  the  sun  (angle  through  which  the  earth  must  turn  to  bring 
the  meridian  of  the  observation  point  directly  under  the  sun) 

6  varies  between  23°  27'  at  the  summer  solstice  (22  June)  and 
-23°  27'  at  the  winter  solstice  (22  December)  and  is  dependent  upon  geographical 
location.  The  hour  angle  varies  by  15°  for  each  hour  difference  from  true 
solar  noon,  at  which  time  it  is  zero.  Thus,  the  insolation  parameter  *  can 
be  determined  with  the  knowledge  of  one  meteorological  observation, 
cloudiness. 

Comparisons  of  the  insolation  parameter  with  measured  values  of 
instantaneous  solar  radiation  are  shown  in  Figure  2.  Each  dot  represents 
an  occurrence.  The  solid  line  represents  the  linear  regression  line  for  this 
comparison.  The  dashed  lines  represent  the  95%  confidence  interval;  i.e., 
any  data  point  chosen  at  random  has  a  95%  probability  of  lying  between  these 
lines.  The  coefficient  of  determination,  r2 ,  for  a  set  of  data  represents  a 
measure  of  the  degree  of  relationship  between  variables.  If  there  were  a 
perfect  correlation  between  the  data,  r  would  equal  1.0,  and  the  data  would 
all  lie  along  one  line.  Note  that  the  linear  regression  line  does  not  really 
adequately  describe  the  data.  The  r2  value  from  the  data  in  Figure  2  is 
0.696. 

Figure  2  represents  a  linear  regression  of  $  and  solar  radiation 
measurements.  A  linear  regression  is  a  first  degree  polynomial  fit. 

Increasing  the  degree  of  the  polynomial  regression  usually  increases  r2,  the 
coefficient  of  determination.  However,  for  the  data  plotted  in  Figure  2,  r2 
for  a  second  degree  polynomial  fit  had  the  value  .703,  essentially  the  same  as 
the  linear  (first  degree)  regression.  The  data  was  too  scattered  to  properly 
fit  an  equation  that  can  be  used  for  prediction  purposes,  but  the  r2  value 
showed  there  Is  a  strong  relationship  between  *  and  this  measurement  of  solar 
radiation. 

To  determine  the  effect  the  albedo  had  on  the  relationship  between  $ 
and  solar  radiation,  the  value  for  A  in  equation  1  was  varied.  Values  between 
0.4  and  0.6  were  chosen.  Although  the  linear  regression  equations  changed, 
no  real  differences  in  correlation  resulted  from  changing  the  albedo  value. 

For  example,  r2  had  the  values  .682  and  .697  for  A  =  0.4  and  A  =  0.6,  respec¬ 
tively. 

To  illustrate  regression  analysis  as  a  prediction  tool,  Figure  3  lists 
measured  radiation  and  the  predicted  radiation  using  equation  4,  the  linear 
regression  equation  derived  from  Figure  2. 


where  <t>  is  defined  in  equation  2. 

Figure  3  shows  the  comparisons  of  equation  4  with  the  instantaneous 
solar  radiation  measurement  taken  at  same  time  of  day  the  cloudiness  was 
estimated  on  7  July  and  21  December.  Figure  3  shows  that  equation  4  fits 
well  with  the  21  December  data,  but  not  the  7  July  data.  The  general  shape 
of  the  7  July  prediction  curve  is  correct  (the  line  predicts  an  increase  in 
solar  radiation  until  12  noon  to  1  p.m,  then  falls  off),  but  the  magnitude 
of  the  estimate  does  not  correspond  to  the  measured  value.  More  than 
anything  else,  this  may  point  out  the  problem  with  comparing  the  measurement 
of  solar  radiation  at  a  point  with  an  insolation  parameter  ( <t> )  intended 
for  a  sufficiently  large  area.  Additional ly ,  comparing  the  instantaneous 
(5-minute  average)  measurement  of  solar  radiation  is  also  risky.  For  example, 
even  though  much  of  the  sky  is  obscured  by  cloud,  the  pyranometer  may  be 
"seeing"  radiation  through  small  openings,  which  would  sense  more  radiation 
than  the  area  in  this  vicinity.  A  time-average  value  of  the  solar  radiation 
should  provide  more  accurate  estimates. 

Figure  4  shows  the  correlation  between  the  insolation  parameter  *  and 
the  average  radiation  over  the  next  60  minutes,  i.e.  ,  how  well  ®  can  be  used 
to  predict  the  upcoming  solar  radiation.  Note  that  the  95%  confidence 
interval  is  narrower  than  that  in  Figure  2,  indicating  a  better  fit.  In  fact, 
the  r2  value  associated  with  Figure  4  is  .792.  This  comparison  would 
indicate  that  a  time-average  of  the  solar  radiation  is  more  reasonable  than 
an  instantaneous  measurement. 

Figure  5  has  been  included  to  show  the  effects  of  changing  the  degree 
of  the  polynomial  for  the  regression  analysis.  Note  the  change  in  slope  of  the 
regression  line  at  the  bottom  left  of  the  graph,  the  area  where  the  observations 
are  most  cluttered.  The  r2  value  for  Figure  5  is  .798,  and  the  regression 
equation  (solid  line  on  Figure  5)  can  be  described  by  equation  5. 

Solar  ( ly/min)  =  .  1177  -  .  9282  *  +  7.486  *2  -  8.461  $3  +  3.194  (5) 

radiation 

Changing  the  albedo  value  in  equation  1  again  had  little  effect  on  the  correlation 
of  these  observations. 

The  insolation  parameter  4>  is  essentially  an  observation  at  a  single 
moment  in  time,  v  will  change  even  in  short  time  spans  because  the  solar 
elevation  constantly  changes.  In  addition,  the  cloud  cover  can  often  change 
rapidly.  To  determine  these  possible  effects,  *  could  be  averaged  over  short 
time  periods  such  as  1  hour.  For  these  purposes  $p  will  be  the  insolation 
parameter  averaged  over  a  period  of  1  hour  previous  to  the  solar  radiation 
measurement. 

For  example,  Figure  6  shows  the  relationship  between  4>p  and  the  instan¬ 
taneous  measurement  of  insolation.  As  noted  earlier,  the  instantaneous 
measurement  of  radiation  is  not  always  a  good  indicator  of  the  amount  of 
radiation  the  general  area  is  receiving.  Note  that  there  is  a  correlation 
between  the  two  measurements ,  but  you  would  not  feel  very  confident  predicting 
the  instantaneous  radiation  using  4>p. 
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Using  ♦p  to  predict  the  average  radiation  should  provide  more 
accurate  estimates.  Figure  7  shows  the  relationship  between  the  insolation 
parameter  averaged  over  the  same  time  period  as  the  solar  radiation.  The 
relationship  is  more  accurate  than  that  evident  from  the  measurement  of  the 
instantaneous  radiation.  The  correlation  of  coefficient  for  the  data  used 
in  Figure  7  is  .772. 

Once  again  employing  *p  as  a  predictor.  Figure  8  shows  the  effects 
of  in  estimating  the  time-lag  solar  radiation.  Here  *p  is  the  average 
insolation  parameter  over  a  given  nour,  and  the  radiation  measurement  is 
the  average  value  for  the  next  hour.  This  relationship  provided  the  best 
correlation,  as  determined  by  the  correlation  coefficient  r? ,  of  any  variables 
analyzed  for  this  study.  The  corresponding  r?  value  for  Figure  8  is  .800. 

It  is  interesting  to  note  that  the  r^  for  this  same  set  of  observations  is 
.813  when  the  albedo  is  defined  to  be  0.6.  Figure  8  shows  that  the  insolation 
parameter  seems  to  be  more  accurate  when  it  is  averaged  over  a  short  time 
period  (~  1  hour)  and  used  as  a  predictor  instead  of  a  measure  of  the  current 
conditions. 

Comparisons  of  the  hourly  delay  insolation  with  two  insolation 
parameters,  *  and  *p,  are  shown  in  Figure  9.  The  measured  radiation  measure¬ 
ments  are  the  same  as  in  Figure  3.  For  the  7  July  case,  *p  provides  the 
better  descriptor  of  the  actual  radiation,  but  it  overestimates  the  radiation 
until  about  1  p.m,  after  which  time  it  underestimates  it.  However,  *  proved 
to  be  the  better  indicator  of  solar  radiation  in  the  21  December  case.  The 
*  curve  tol lows  the  measured  radiation  curve  closely.  Note  that  for  the 
summer  case  both  predictors  overestimate  the  measured  radiation  for  the 
first  half  of  the  day  and  underestimate  it  for  the  second  half,  and  the 
reverse  is  true  for  the  winter  case.  There  is  no  apparent  reason  for  this 
pattern. 

Table  10  summarizes  the  different  ways  in  which  the  insolation 
parameter  and  solar  radiation  were  compared,  along  with  the  associated 
correlation  coefficient. 

Table  10.  Summary  of  Comparisons  between  the  Insolation  Parameter  and 
Measured  Solar  Radiation 


Insolation  parameter  Albedo  Solar  radiation  r?  (correlation  coefficient) 


*  0.5  Instantaneous  .696 
*p  0.5  Instantaneous  .724 
4>p  0.5  Average  hourly  .7/2 

*  0.5  Hourly  delay  .792 
4>p  0.5  Hourly  delay  .8UU 
$n  0.6  Hourly  delay  .813 


S 


There  are  two  possible  reasons  that  the  correlations  summarized 
in  Table  10  are  not  higher.  First,  all  solar  radiation  measurements  were 
made  with  one  stationary  pyranometer,  essentially  a  fixed  point.  The 
insolation  parameter  4>  is  intended  to  estimate  the  strength  of  solar  radiation 
for  a  general  area.  The  problem  of  cloudiness  over  a  general  area  and  that 
at  a  point  were  discussed  previously.  There  is  no  doubt  that  averaging  at 
least  two  measurements  of  solar  radiation,  adequately  spaced,  would  improve 
this  relationship.  Secondly,  K  is  not  really  constant  as  assumed  in  equation 
1.  K  is  a  function  of  the  solar  constant,  which  is  quite  constant,  and 
atmospheric  transmission,  which  fluctuates  with  time.  The  amount  of  scattering 
and  absorption  varies  with  sky  conditions,  including  synoptic  (large-scale) 
weather  patterns  and  atmospheric  pollutants. 

Another  way  to  determine  the  relationship  between  two  parameters 
is  by  using  the  joint  occurrence  table,  showing  the  frequency  of  occurrence 
of  one  parameter  with  another.  Table  11,  which  categorizes  the  same  obser¬ 
vations  used  in  Figure  2,  is  an  example.  Table  11  details  how  the  sum  of 
all  occurrences  of  one  parameter  corresponds  with  the  other.  For  example, 

18.7%  of  all  observations  occurred  with  .30  <  .49  and  .30  <  SR  <  .69. 

Also,  when  .30  <  *  <  .49,  65.4%  of  the  solar  radiation  measurements  fell 

AS.7\ 

between  .  30  and  .  69  [ -  1=  .  654  . 
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Table  11.  Joint  Occurence  of  Insolation  Parameter  $and  Instantaneous 
Solar  Radiation 

(214  observations,  A  =  0. 5) 


<  .3 

Instantaneous  solar  radiation  (ly/min) 
.30-69  .70-1.09  >  1.1 

<.3 

29.9 

5.6 

0 

0 

35.5 

.30-.  49 

6.1 

18.7 

3.3 

0.5 

28.6 

.50-.  69 

1.4 

4.7 

9.8 

3.3 

19.2 

>.70 

0 

1.8 

5.6 

9.3 

16.7 

37.4 

18. 1 

nn 

100.0 

Tables  12  through  15,  corresponding  to  Figures  4,  6,  7,  and  8, 
respectively,  are  shown. 


Table  12.  Joint  Occurrence  of  Insolation  Parameter  $  and  Average 
Hourly  Delay  Radiation  (187  observations,  A  =  0.5) 


E 


*p 

<.3 

Average  hourly  delay  radiation  i 
.30-. 69  .70-. 99 

ly/min) 

>1.0 

<.3 

27.3 

5.3 

0 

0 

.30-. 49 

7.5 

20.9 

2.7 

0 

.50-. 69 

0.5 

4.3 

9.1 

3.2 

>  .70 

0 

1.1 

3.7 

14.4 

35.3 

31.6 

15.5 

17.6 

9.2 


100.0 


Table  13.  Joint  Occurrence  of  Insolation  Parameter  4>p  and  Instantaneous 
Solar  Radiation  (187  observations,  A  =  0.5) 


Instantaneous 

solar  radiation 

Ty/min) 

— 

<.3 

.30-. 69 

.70-1.09 

>1.1 

<.3 

26.2 

7.5 

0 

0 

.30-. 49 

5.3 

16.6 

2.7 

0 

.50-. 69 

1.6 

5.3 

13.9 

4.8 

>  .70 

0 

1.1 

4.8 

10.2 

33.1 

30.5 

21.4 

15.0 

100.0 


Table  14.  Joint  Occurrence  of  Insolation  Parameter  $p  and  Average  Hourly 
Radiation  (187  observations,  A  =  0.5) 
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Table  15.  Joint  (a  .  rrence  of  Insolation  Parameter  fcp  and  Average  Hourly 
Delay  Radiation  (187  observations,  A  =  0.5) 


Average  hourly  delay  radiation  (ly/min) 
.30-. 69  .70-1.09  >1.1 


<.3 

27.3 

6.0 

0.5 

0 

33.8 

.30-. 49 

7.0 

16. 0 

1.6 

0 

24.6 

.50-. 69 

0.5 

4.8 

16.0 

4.3 

25 . 6 

>  .70 

0 

0.5 

3.2 

12.3 

16.0 

34.8 

27.3 

21.3 

16.6 

100.0 

Had  all  the  observations  been  located  along  the  diagonal  running 
from  top  left  to  bottom  right  in  each  table,  this  method  of  comparison 
would  have  implied  a  "perfect"  relationship  between  the  variables.  The 
sum  of  the  observations  along  the  diagonals  ranged  from  66.9%  (Table  13)  to 
71.6%  (Table  15)  and  71.7%  (Table  12).  The  measurement  of  insolation  in 
the  form  of  hourly  averages  was  found  to  have  the  highest  percentage  of 
observations  along  the  diagonals,  while  the  instantaneous  measures  provided 
the  smallest  percentage  of  observations  along  the  diagonals. 

Similar  to  the  results  from  the  regression  analyses  discussed  earlier, 
the  albedo  value  in  Tables  11-15  had  little  effect  on  the  relationship  between 
the  insolation  parameter  $  and  the  insolation  measurements. 

9.  DATA  COMPARISON  -  SMITH  ALGORITHM 

For  comparative  purposes,  a  previous  study  of  solar  radiation  has 
been  included.  Smith  developed  an  algorithm  (Figure  10)  for  determining 
the  Pasquill  Stability  Parameter  from  the  direct  measurement  of  solar  radiation 
and  wind. 5  Figure  10  includes  the  stability  determination  for  the  daytime 
only.  Smith  included  a  provision  for  nighttime  stability  as  determined  by 
the  upward  heat  flux  at  the  surface,  but  this  has  not  been  included  since 
nighttime  measurements  are  not  used  in  this  study. 

The  most  noticeable  feature  of  Figure  10  is  the  determination  of 
stability  on  a  continous  scale.  Recall  that  the  P-T  system  categorized  the 
daytime  stability  parameter  into  one  of  four  distinct  values  ranging  from  A  to 
D.  The  Smith  Algorithm  allows  split  categories.  For  example,  the  stability 
parameter  2.2  is  towards  the  more  unstable  end  of  category  C.  (These  split 
categories  are  important  in  determining  az  (x),  the  vertical  standard 
deviation  of  the  chemical  cloud,  a  measure  of  the  vertical  growth.  This 
parameter  determines  the  downwind  concentration  or  dosage  of  the  chemical 
cloud.) 

Two  patterns  are  evident  in  Figure  10.  As  the  solar  radiation 
increases,  assuming  a  constant  wind  speed,  the  stability  parameter  becomes 
more  unstable.  Secondly,  keeping  solar  radiation  constant  and  increasing 
the  wind  speed  pushes  the  stability  parameter  toward  neutral  ,  the  most  stable 
condition  possible  during  the  daytime.  These  same  patterns  existed  in  the 
P-T  stability  algorithm  (Table  3). 


Windspeed (at  10m) 


in  tj* 


'  i 

i  i 

*  i 

i  i 

•  i 

i  ! 


i  ;  i  > 
I  1 

i  1  /  : 

I  :  I  : 

i  1  ; 

I  1 

/  :  ! 


!  i  i  .  // 


;  /  #  '  /  / 
;  /  /  '  /  / 

i  i  1  '  i  •  j 

lit  '  /  /  / 

i  i  '  ■  /;  / 

>  I  '  /  /  / 
l  ;  '  • // / 

i  / ,  //  / 


SOLAR  RADIATION  CLY/MIN) 


uu IW  W  VV  V V n  ■-'  M  ■* ^  V -.v-v-y >  V * v»  v«  VTV * ;'V'v«VTv»\'vin.1vwww\'n 


Note  that  the  neutral  stability  is  assigned  the  value  3.6  in  Figure  10. 
If  numerical  values  were  to  be  assigned  to  the  P-T  categories,  A  would  equal 
1.0,  B  =  2.0,  C  =  3.0,  and  0  (neutral)  a  4.0.  The  corresponding  values  from 
Figure  10  are  0.6,  1.6,  2.6,  and  3.6,  respectively,  which  represent  the  center 
of  each  P-T  stability  category.  Thus,  to  compare  stability  parameters,  add 
0.4  to  the  Smith  value  to  obtain  the  P-T  stability  parameter. 

Smith  also  Included  a  provision  for  determining  stability  based 
on  the  upward  heat  flux  (UHF),  the  transfer  of  heat  from  the  earth's 
surface  (see  Section  2).  The  UHF  could  be  measured  directly  or  it  could 
be  estimated  using  equation  5  and,  if  measured  directly,  could  be  used  for  a 
nighttime  stability  indicator. 5 


UHF  =  0.4  (SR  -  .1433) 


(6) 


where,  UHF  and  solar  radiation  have  units  (ly/min), 


Assuming  the  measurement  of  UHF,  equation  5  allows  the  calculation 
of  stability  on  a  continuous,  daily  scale.  Neither  the  UHF  scale  nor  the 
nighttime  stability  curves  are  shown  in  Figure  10. 


Table 

solar  radiation 
discrete  stabil 
on  a  continuous 
not  recorded, 
average.  Table 
corresponding  P 
height,  and  sol 


16  lists  the  P-T  stability  category  as  determined  using  the 
data  summarized  in  this  report.  Table  16  considers  only 
ity  parameters,  since  the  stability  parameter  was  not  recorded 
scale.  The  table  Includes  interpolations  where  data  was 
The  solar  radiation  measurement  represents  a  1-hour  time 
16,  then,  summarizes  the  average  hourly  radiation  with  the 
-T  stability  Index  as  determined  by  the  cloud  cover,  cloud 
ar  elevation. 


Table  16.  Pasquill  Stability  Category  Based  on  the  Direct  Measurement  of 
Solar  Radiation 


Table  16  lists  the  most  frequent  occurrence  of  the  stability  index 
for  the  corresponding  wind  speed  and  solar  radiation  measurements.  Some 
split  categories  existed  and  are  denoted  by  a  (-). 

The  same  general  trends  shown  in  Figure  10  are  seen  in  Table  16: 

(a)  increasing  the  solar  radiation  as  the  wind  speed  remains  constant 
decreases  the  stability  parameter  (more  unstable)  and  (b)  increasing  the 
wind  speed  as  the  solar  radiation  remains  constant  increases  the  stability 
parameter  (toward  neutral). 

The  greatest  range  of  stability  occurred  when  the  anemometer  recorded 
calm  conditions.  The  stability  parameter  varied  from  D  to  A  over  a  range  of 
only  0.6  ly/min.  It  Is  Interesting  to  note  that  the  occurrence  of  stability 
A  was  virtually  confined  to  calm  winds.  Outside  the  zero  wind  category, 
stability  A  appears  on  Table  16  only  when  the  wind  speed  equaled  1  m/sec  and 
solar  radiation  measured  at  least  1.2  ly/min.  Note  that  repetition  allowed 
the  combination  of  some  radiation  categories,  i.e.,  .7-. 9  and  >  1.2  ly/min. 

In  general,  the  data  in  Table  16  compares  favorably  with  that  in 
Figure  10.  For  calm  winds,  the  stability  parameters  match  except  for  the 
0.4  and  0.6  Insolation  categories,  where  In  Table  16  one  stability  category 
closer  to  the  unstable  end  Is  predicted.  For  a  2-m/sec  wind  speed,  the  curves 
agree  quite  well  except  In  the  range  .85-1.00  ly/min,  where  in  Figure  10 
an  A  stability  Is  predicted  while  In  Table  16  B  stability  Is  indicated. 

At  the  3-m/sec  wind,  the  curves  agree  up  to  0.5  ly/mln  insolation,  but 
beyond  this  value  (Figure  10)  a  more  unstable  category  is  predicted. 

Finally,  In  Table  16  a  neutral  stability  category  up  to  at  least  0.4  ly/min 
Is  predicted,  while  In  Figure  10  a  neutral  category  only  to  approximately  0.2 
ly/min  Is  indicated.  Due  to  Insufficient  data  involving  higher  wind  speeds, 
comparisons  between  the  two  systems  could  not  be  made.  However,  for  all 
observations  where  comparisons  could  be  made,  the  two  systems  never  differed 
by  more  than  one  stability  category. 

One  point  concerning  this  comparison  should  be  emphasized.  The  data 
used  to  generate  Table  16  was  collected  In  the  vicinity  of  trees  and  buildings 
near  an  Inlet  of  the  Chesepeake  Bay.  The  wind  barriers  and  large  water  body 
affect  both  the  local  circulation  of  winds  and  the  surface  heating  (not  how 
much  Insolation  received  per  se,  but  how  much  retained),  the  two  governors 
of  atmospheric  stability.  On  the  other  hand,  Figure  10  applies  to  open 
terrain.  Therefore,  Table  16  may  have  to  be  revised  slightly  if  it  is  to 
apply  to  a  generalized,  mid-latitude  temperate  climate. 

Table  17  has  been  included  for  comparative  purposes,  displaying 
the  Pas  quill  stability  category  as  a  function  of  the  insolation  parameter 
•p  and  wind  speed.  Table  17  was  developed  from  the  same  observations  as 
Table  16,  and  *p  represents  a  1-hour  average  for  the  same  observations.  A 
direct  comparison  to  Figure  10  Is  not  possible,  but  the  same  trends  of 
stability  found  In  Figure  8  and  Table  3  also  exist  here. 


Table  17.  Pasqulll  Stability  Parameter  Based  on  ap 
•p  (dimensionless) 


wind  speed 


(m/sec) 

0-0.2 

.3 

.4 

.5 

.6 

.7 

>  .80 

0 

D 

C 

B 

A 

A 

A 

A 

1 

D 

C 

B 

B 

B 

A 

A 

2 

D 

C 

C 

C 

B 

B 

A 

3 

D 

C 

C 

C 

C 

B 

4 

D 

D 

C 

5 

0 

D 

>  6 

D 

D 

10.  CONCLUSIONS 

Two  subjective  estimates  of  solar  radiation  have  been  compared 
with  pyranometer  measurements.  Both  the  P-T  and  Insolation  parameter  *  methods 
displayed  a  reasonably  good  fit  with  the  field  data.  In  comparing  the  two 
subjective  estimates,  the  insolation  parameter  *  was  judged  to  perform 
better  than  the  P-T  system. 

The  relationship  between  all  variables  improved  when  time-averaged 
quantities  were  considered. 

Three  factors  negatively  Influenced  the  performance  of  the  Insola¬ 
tion  parameter  *: 

•  The  proportionality  factor  dependent  on  atmospheric  transmission 
Is  not  constant.  The  amount  of  solar  radiation  allowed  through  the  atmosphere 
can  vary  greatly  over  the  range  of  synoptic  (large-scale)  weather  conditions. 

•  The  average  cloud  albedo  used  does  not  apply  to  all  cloud  types. 

•  The  solar  radiation  measurements  were  collected  at  a  single  point. 
Fluctuations  in  Insolation  due  to  Irregular  cloud  distributions  can  be 
expected  at  a  point,  which  may  not  accurately  describe  the  solar  radiation 
being  received  over  the  general  area.  It  Is  reasonable  to  expect  these 
fluctuations  to  be  smoothed  through  the  averaging  of  Insolation  at  several, 
well -spaced  points. 

Several  factors  may  have  negatively  Influenced  the  performance  of 
the  P-T  system: 

•  Cloud  cover  and  cloud  height  measurements  were  made  subjectively. 
Cloud  height  was  estimated  based  on  the  Identity  of  cloud  type.  However, 
several  cloud  types  are  similar  In  appearance  and  difficult  to  distinguish. 

In  addition,  such  atmospheric  conditions  as  haze  can  make  It  difficult  to 
quantify  cloud  cover  and  height. 

•  There  Is  reason  to  believe  that  the  anemometer  recordings  of 
wind  speed  are  low,  l.e.,  the  actual  wind  speeds  were  higher  than  that  recorded. 
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•  The  solar  radiation  measurements  were  collected  at  a  single 
point,  while  the  P-T  system  defines  atmospheric  stability  for  a  generalized 
area.  There  is  no  doubt  this  accounts  for  some  error  in  comparison  of  the 
two  quantities. 

Algorithms  have  been  provided  that  select  the  Pasquill  Stability 
Category  based  on  the  direct  measurement  of  solar  radiation.  These  algorithms 
eliminate  the  need  to  subjectively  estimate  insolation.  Where  the  direct 
measure  of  solar  radiation  is  not  possible,  the  insolation  parameter  $  and 
P-T  systems  have  been  shown  to  be  adequate  methods  for  characterizing  solar 
radiation. 

The  ultimate  use  of  all  atmospheric  stability  predictors  is  to 
estimate  the  dispersion  rates  of  vapors,  aerosols,  and/or  smokes  under 
those  meteorological  conditions.  It  is  recognized  that  this  study  does  not 
predict  the  behavior  of  a  chemical  cloud  introduced  to  the  atmosphere,  but 
it  provides  a  slightly  simpler  method  to  determine  atmospheric  stability. 

11.  RECOMMENDATIONS 

Data  from  time  periods  other  than  those  recorded  here  should  be 
collected  and  analyzed  to  determine  if  the  results  presented  in  this  study 
contrast  significantly.  It  is  expected  that  the  average  solar  radiation 
measurements  associated  with  stability  category  D  will  prove  to  be  higher 
when  more  data  has  been  collected.  In  this  study,  most  incidences  of  stability 
D  occurred  near  sunrise  or  sunset,  when  insolation  had  to  be  low.  However, 
stability  D  can  also  occur  during  miday  under  heavy  cloud  cover.  The  solar 
radiation  reaching  the  surface  should  be  higher  under  these  circumstances. 

Data  should  be  collected  using  at  least  two,  well-spaced  pyranometers. 

Data  should  be  reduced  automatically.  A  significant  portion  of 
time  tor  this  study  was  spent  reducing  data. 

Measurements  of  surface  heat  flux  could  be  used  to  determine  its 
relationship  to  surface  insolation.  These  measurements  could  then  be  used  to 
determine  atmospheric  stability  during  the  daytime  or  nighttime.  Solar 
radiation  measurements,  obviously,  can  only  be  used  to  determine  daytime 
stability. 
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